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Theoretical Characterization of the Reaction 

CH S + OH -> CH 3 OH -> products : 

The X CH 2 + H 2 0, H 2 + HCOH, and H 2 + H 2 CO Channels 

Stephen P. Walch“ 

ELORET Institute 
Palo Alto, Ca. 94303 

Abstract. The potential energy surface (PES) for the CH 3 OH system has been char- 
acterized for the X CH 2 + H 2 0, H 2 + HCOH, and H 2 + H 2 CO product channels 
using complete-active-space self-consistent-field (CASSCF) gradient calculations to 
determine the stationary point geometries and frequencies followed by CASSCF/ 
internally contracted configuration-interaction (CCI) calculations to refine the ener- 
getics. The x CH 2 + H 2 0 channel is found to have no barrier. The long range inter- 
action is dominated by the dipole-dipole term, which orients the respective dipole 
moments parallel to each other but pointing in opposite directions. At shorter sep- 
arations there is a dative bond structure in which a water lone pair donates into the 
empty a" orbital of CH 2 . Subsequent insertion of CH 2 into an OH bond of water 
involves a non-least-motion pathway. The H 2 + HCOH, and H 2 + H 2 CO pathways 
have barriers located at -5.2 kcal/mol and 1.7 kcal/mol, respectively, with respect 
to CH 3 + OH. From comparison of the computed energetics of the reactants and 
products to known thermochemical data it is estimated that the computed PES is 
accurate to i 2 kcal/mol. 

“Mailing Address: NASA Ames Research Center, Moffett Field, CA 94035. 
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I. Introduction 

The CH 3 + OH reaction has at least six possible product channels: 


CH 3 + OH + M - CH 3 OH + M (la) 

CH 3 + OH — ► CfTaOJ? + H ( U ) 

CF 3 + OH -> CH 3 0 + H ( lc ) 

CHz + OH ^ CH 2 + H 2 0 0 - d ) 

CH 3 + OH — > HCOH + H 2 ( le ) 

CH 3 + OH H 2 CO + H 2 (!/) 


The role of reaction Id has been controversial. The room temperature rate 
for the reverse of reaction Id has been measured by Hatakeyama et al. [1] 
as ~ 3 x 1CT 12 cm 3 molecule -1 cm -1 , while Hack et al. [2] obtained 
3.5 x 10 -n cm 3 molecule -1 cm -1 , or about an order of magnitude faster. The latter 
rate is approximately gas kinetic and implies no barrier. Using currently accepted 
heats of formation [3], and a singlet-triplet splitting for methylene of 9 kcal/mol, 
reaction Id is exothermic by 0.7 kcal/mol at OK. Thus, the forward reaction is also 
expected to be very fast. 

Dean and Westmoreland [4] have used a variant of Rice-Ramsperger-Kassel- 
Marcus (RRKM) theory called QRRK theory to model the product distributions 
in the CH 3 + OH reaction. In this work the parameters for lb and lc were based 
on estimated rates for the reverse reaction. The parameters for reaction If were 
taken from calculations [5] and the rate for reaction Id was taken from Ref. 1. 
According to these studies, at room temperature and moderate pressure CH 3 OH is 
the dominant product, while at flame temperatures reaction lb is thought to take 


2 


over. The HCOH + H 2 channel does not appear to have been considered in these 
studies, though theory [5] indicates essentially no barrier with respect to CH 3 + 
OH. This study indicated that production of 1 CH 2 (reaction Id) is a minor channel. 
By contrast a model proposed by Pilling and coworkers [6] which makes use of the 
rate for the reverse of reaction Id due to Hack et al. [2] indicates that reaction Id 
is the dominant channel above room temperature. 

Recently Smith [7] has also reported RRKM calculations for CH 3 + OH. These 
calculations as well as the work of Pilling et al. [6] have indicated a need for more 
accurate potential energy surface (PES) information for the 1 CH 2 + H 2 0, H 2 + 
HCOH, and H 2 + H 2 CO product channels. The most accurate previous theoretical 
study of these channels in the CH 3 OH system was carried out by Harding, Schlegal, 
Krishnan, and Pople [5] using Mpller-Plesset perturbation theory with a 6-311G** 
basis set. Although these calculations were carefully carried out, there is probably 
considerable uncertainty in the energetics; by current standards, both the basis set 
and treatment of electron correlation can be improved upon. More recently the bond 
dissociation energies of CH 3 OH have been computed by Bauschlicher, Langhoff, and 
Walch [8] using the modified coupled-pair functional method. Similar calculations 
were carried out by Pople and co-workers [9] using the G2 method, which includes 
some empirical corrections. These calculations accurately determined the heats of 
formation of the CH 3 0 and CH 2 OH species, but they did not consider the portions 
of the PES leading to the 1 CH 2 + H 2 0, H 2 + HCOH, and H 2 + H 2 CO product 
channels. Thus, these regions of the PES are reexamined here. 

Qualitative features of the potential energy surfaces are discussed in Sec. II, the 
computational method is discussed in Sec. Ill, the results are presented in Sec. IV, 
and the conclusions are given in Sec. V. 
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II. Qualitative Features. 

. , . . , _ 3d gr ound state which will be drawn as 

CH 2 has two low-lying states, lne 1 g 



has two orthogonal high- spin 


coupled singly occupied orbitals, while the 


^ A. 1 S t £lt 0 , 


which is « 9 kcal/mol higher, is drawn as 



he., the two orbitals corresponding to the C lone pair are singlet paired. If these 
orbitals are solved for self consistently in a generalized valence bond wavefunc 
[101 the overlap integral between them is * 0.7. Thus, the singlet state of CH, may 
be characterized as a singlet biradicU, but the orbitals of the lone pair have a sub- 
stantial overlap, which must be maintained while inserting into a single bond if the 
process is to occur without a substantial barrier. In terms of multiconfigurational 
self-consistent-field (MCSCF) theory, the biradical character in CH 2 arises because 
of a near degeneracy effect between an sp hybrid lone pair and an empty 2p-Hke a" 


orbital. 

The insertion 
This process is 


of >CH 2 into a bond pair occurs via a non-least-motion pathway, 
completely analogous to the non-least-motion addition of the ’II 
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state of CH to H 2 , as discussed by Dunning and Harding [10]. In that case the CH 
and H 2 approach each other with the bond axes parallel and for this orientation 
the C lone pair and H 2 bond pair orbitals are able to evolve into two CH bonding 
orbitals without breaking either bond. These orbital changes are consistent with 
the orbital phase continuity principle arguments made by Goddard [11]. 

In the case of ^Ho + H 2 0 the long range interaction is dominated by the dipole- 
dipole term, which leads to an initial approach with the dipoles parallel to each other 
but pointed in the opposite direction. At shorter separation there is a dative bonded 
structure where one lone pair on H 2 0 donates into the empty C 2p-like orbital of 
CH 2 . This structure is a minimum on the PES. Insertion of CH 2 into one of the 
H 2 0 bonds requires rotating H 2 0 such that one OH bond is in the plane defined by 
O and the bisector of L HCH. This orientation is similar to the orientation favored 

for non-least-motion insertion of CH into H 2 • 

Insertion of 1 CH 2 into H 2 0 involves very little barrier. The insertion of hydroxy 
methylene (HCOH) into H 2 , however, is found to exhibit a barrier of 13.9 kcal/mol. 
The lower reactivity of HCOH as compared to CH 2 results from delocalization of 
an O lone pair into the empty C 2p orbital. The resultant exclusion effect reduces 
the 2s — > 2p near-degeneracy effect and thus results in less biradical character in 

the substituted carbene. 

III. Computational Details. 

Two different basis sets were used in this work. For the CASSCF gradient calcu- 
lations the polarized double-zeta set of Dunning and Hay [12] was used. The basis 
set for C and 0 is a (9s5p)/[3s2p] basis augmented by a single set of 3d functions 
with exponents of 0.75 and 0.85 for C and 0, respectively. The H basis is (4s)/[2s] 
augmented with a single set of 2p functions with exponent 1.00. The basis set 
used in the Cl calculations is the Dunning correlation consistent tnple-zeta double- 
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polarization basis set [13]. This basis is [4s3p2dlf] for C and O and [ 3 s 2 pld] for H 
and is described in detail in Ref. 13. 

In the CASSCF calculations the electrons in the bonds which are being made 
or broken were included in the active space. In the case of the 1 CH 2 + H 2 0 and 
H 2 CO + H 2 channels, there were six active electrons and six active orbitals. The 
correlated electrons correspond to the CO bond, the OH bond, and one CH bond 
of CH 3 OH and evolve to the C lone pair of CH 2 and the two OH bond pairs m the 
case of j CH 2 + H 2 0 and the CO a and tt bonds and the H 2 bond in the case of 
H 2 CO + H 2 . The remaining 12 electrons are inactive in the CASSCF calculation. 
In the case of the HCOH + H 2 channel, there were four active electrons and four 
active orbitals. The correlated electrons correspond to two CH bonds m CH 3 OH 
and to the C lone pair in HCOH and the H 2 bond pair in the HCOH + H 2 limit. 
The remaining 14 electrons were inactive in the CASSCF calculation. In generating 
the set of reference configurations, no more than two electrons were permitted m 
the weakly occupied CASSCF orbitals. All but the 0 Is and C Is electrons were 
correlated in the CCI calculations. 

The CASSCF/gradient calculations used the SIRIUS/ABACUS system of pro- 
grams [14], while the CCI calculations were carried out with MOLPRO [15,16]. 
Most of the calculations were carried out on the NASA Ames Cray Y-MP; although 
some of the CCI calculations were carried out on the NAS facility Y-MP. 

IV. Discussion. 

Table I shows computed energetics for the portions of the CH 3 OH surface which 
were considered in this work. The computed energetics were obtained from the CCI 
energies (including the multireference analogue of the Davidson correction [17]) plus 
the zero-point energies from the CASSCF calculations. (See Tables Ila-IIc.) Thus, 
these energetics should be compared to experimental values at 0 K. The computed 
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frequencies and rotational constants for the saddle points for CH 2 0 + H 2 , HCOH 
+ H 2 , and 1 CH 2 + H 2 0, ( denoted as CH 2 0-H 2 , HCOH-H 2 , and CH 2 -H 2 0, re- 
spectively) and for the CH 2 + H 2 0 dative bonded structure (denoted as CH 2 .H 2 0) 
are also given in Table III. The stationary point corresponding to CH 2 .H 2 0 is a 
minimum on the PES, but there is one very small frequency ( 26 cm : ), which 
corresponds to a hindered rotation of the CH 2 and H 2 0 with respect to each other. 

The computed energy separations discussed here, in each case, involve breaking 
two bonds and forming two new bonds; thus, the errors in the individual bond 
strengths cancel and the computed energetics are expected to be accurate. However, 
in Ref. 2 it was shown that, for calculations of about the same quality as reported 
here, the error in the C-0 bond strength in CH 3 OH is 6.5 kcal/mol. Thus, in order 
to compute energies with respect to CH 3 + OH, the experimental 0 K value of 
90.2 kcal/mol [18] was used for the C-0 bond strength. The locations of the H + 
CH3O and H + CH 2 OH asymptotes were taken as the best-estimate values from Ref. 
2. This places H + CH 3 0 and H + CH 2 OH at 14.8 kcal/mol and 6.0 kcal/mol above 
CH 3 + OH, respectively. The experimental locations of 1 CE 2 + H 2 0 and CH 2 0 + 
H 2 with respect to CH 3 OH were derived from the JANAF [3] heats of formation of 
H 2 , CH 2 0, CH 3 , OH, H 2 0(g), and 3 CH 2 , plus a singlet-triplet splitting in CH 2 of 
9.0 kcal/mol and the value for the C-0 bond strength in CH3OH given above. 

From Table I it is seen that the computed CH3OH -* 3 CH 2 + H 2 0 separation 
is 0.7 kcal/mol smaller than experiment, while the computed CH 3 OH -> CH 2 0 + 
H 2 separation is 0.6 kcal/mol smaller than experiment. It is also seen that the 
computed results of Harding et al. [5] are 5.4 kcal/mol larger and 2.3 kcal/mol 
smaller, respectively, for the same separations. The computed barrier heights for 
CH 3 OH - HCOH + H 2 and CH 3 OH - CH 2 0 + H 2 obtained in Ref. 5 are 
6.0 kcal/mol and 4.6 kcal/mol larger than the barrier heights obtained in the present 
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work. These differences presumably reflect a combination of the larger basis set 
and more extensive correlation treatment in the present calculations. Based on the 
comparison to experiment for the 1 CH2 + H2O and CH2O + H2 asymptotes it is 
reasonable to assign error bars of ~ ± 2 kcal/mol to the computed energetics in 
this work. 

The computed energetics for the CH3OH system are also shown schematically 
in Fig. 1 . From Fig. 1 it is seen that only the CH2O + H2 channel exhibits a 
barrier ( 1.7 kcal/mol.) with respect to CH 3 + OH. The H 2 + HCOH channel has 
a barrier, but it is below CH3 + OH. The 1 CH2 + H2O channel has no barrier and 
is computed to be 1.4 kcal/mol below CH3 + OH. The CH 3 0 + H and CH 2 OH + 
H channels are endoergic with respect to CH 3 + OH, but no intermediate barriers 
are expected. Thus, from this work the 1 CH 2 + H 2 0 , HCOH + H 2 , and CH 2 0 + 
H 2 channels are all found to be accessible from CH 3 + OH. 

Fig. 2 shows the energy as a function of r C o for the addition of 1 CH 2 to H 2 0 , 
while the computed energies are given in Table IV. Two stationary points have been 
located along this minimum energy path. At long rco the dominant interaction 
is dipole-dipole, which results in an orientation with the dipole moments of the 
approaching molecules parallel to each other but pointing in opposite directions, 
as illustrated in Fig. 2. At shorter r C o there is a minimum on the PES with 
C, symmetry, followed by a saddle point with no symmetry for insertion of 1 CH 2 
into an OH bond of water. In order to characterize the minimum energy path, a 
gradient calculation was carried out starting at the saddle point and proceeding 
toward the minimum. CCI calculations were then carried out at the geometries 
corresponding to each step on the walk. This calculation defines the minimum 
energy path between the saddle point and minimum. It is not possible to do the 
same calculation for the portion of the PES connecting the minimum and the 1 CH 2 
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+ H 2 0 asymptote, since there is no well-defined way to follow the gradient uphill. 
In order to characterize this portion of the PES, calculations were carried out with 
the CH 2 and H 2 0 molecules fixed at their equilibrium geometries and oriented with 
the planes of the molecule parallel to each other and the CO bond perpendicular 
to both molecular planes in the orientation shown in Fig. 2. For this geometric 
orientation, rco was varied and the resulting energies at the CCI level are also 
included in Fig. 2. 

The main features of Fig. 2 are a shallow minimum followed by a small barrier 
to formation of CH 3 OH. However, the barrier is below the X CH 2 + H 2 0 asymp- 
tote and therefore the bottleneck on the vibrationally adiabatic curve is expected 
to occur in the entrance channel region. The main feature responsible for the en- 
trance channel bottleneck is the building in of bending modes, which arise from 
electrostatic (dipole-dipole, dipole-quadrapole, and quadrapole-quadrapole) inter- 
actions. In order to define this interaction, dipole and quadrapole moments (about 
the center of mass) were computed for X CH 2 and H 2 0 and are given in Table V. The 
experimental values [19] for H 2 0 (in a.u.) are p. = 0.73, Q** = 1.955, Q yy = -1.859, 
and Q zz = -0.097; which are in reasonable agreement with the computed values. 

V. Conclusions. 

The potential energy surface (PES) for the CH3OH system has been charac- 
terized for the x CH 2 + H 2 0, H 2 + HCOH, and H 2 + H 2 CO product channels 
using complete-active-space self-consistent field (CASSCF) gradient calculations to 
determine the stationary point geometries and frequencies followed by CASSCF/ 
internally contracted configuration-interaction (CCI) calculations to refine the en- 
ergetics. 

The H 2 + H 2 CO, and H 2 + HCOH pathways have barriers located at 1.7 kcal/mol 
and - 5.2 kcal/mol with respect to CH 3 + OH. The X CH 2 + H 2 0 channel is found to 
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have no barrier in the absence of vibrational zero-point effects. However, the long 
range-interaction is dominated by a dipole-dipole term and the zero-point effects 
due to this interaction are expected to lead to a bottleneck on the vibrationally 
adiabatic minimum energy path. The 1 CH 2 + H 2 O asymptote is computed to be 
1.4 kcal/mol below CH 3 + OH. Thus, all three of these channels are expected to be 
accessible at moderate temperatures. 

From comparison of the computed energetics of the reactants and products to 
known thermochemistry it is estimated that the computed PES is accurate to 
dh 2 kcal/mol. 

Acknowledgement. SPW was supported by NASA cooperative agreement number 
NCC2-478. 
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Table I. Computed Energetics for CH 3 0H a,i> 



calc. 

exp(OK) 

Harding et al. 

CH 3 + OH 


90.2 


HCOH-H 2 

85.0 


91.0 

HCOH + H 2 

71.1 


71.1 

1 cb 2 + h 2 o 

88.8 

89. 5 C 

94.9 

ch 2 o-h 2 

91.9 


96.5 

ch 2 o + h 2 

18.1 

18.7 

16.4 

CH 3 OH 

0.0 

0.0 

0.0 

0 CASSCF/CCI with 

a [4s3p2dlf/3s2pld] basis set. 


b Relative energies in 

kcal/mol (including zero-point 

energy). 

c Using 3 CH 2 -» J CH 

2 separation of 9.0 kcal/mol. 
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Table Ila. Computed energies and zero-point corrections. 


CH 3 OH -+ CH 2 0 + H 2 



Energy 0 

zero-point energy** 

A E' 

CH 30 H 

-115.51397(-.54566) 

0.05258 

0.0 

CH 2 0-H 2 

-115. 35603(-. 39071) 

0.04407 

91.9 

ch 2 o+h 2 

- 115. 4721 5(-. 50189) 

0.03768 

18.1 


. Energy in E„. The first energy is the CCI energy, while the energy in parenthesis 
includes a multi-reference Davidson correction and is with respect to -115. Eh- 
b zero-point energy in Eh- 

c relative energy in kcal/mol including zero-point energy and a multi-reference 
Davidson’s correction. 
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Table lib. Computed energies and zero-point corrections. 



CH30H - 

D CH 2 + H 2 O 



Energy a 

zero-point energy** 

A E c 

CH 3 OH 

-115. 51397(-. 54566) 

0.05258 

0.0 

CH 2 -H 2 0 

-11 5.36931 (-.40561) 

0.04366 

82.3 

ch 2 +h 2 o 

-115.35989(-. 39040) 

0.03884 

88.8 


* Energy in E H - The first energy is the CCI energy, while the energy in parenthesis 
includes a multi-reference Davidson correction and is with respect to -115. Eh- 
b zero-point energy in Eh- 

c relative energy in kcal/mol including zero-point energy and a multi-reference 
Davidson’s correction. 
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Table lie. Computed energies and zero-point corrections. 


CH z OH -» HCOH + H 2 



Energy 0, 

zero-point energy^ 

A E 

CH 3 OH 

-115. 50589(-. 54231) 

0.05258 

0.0 

HCOH-H 2 

-115.35803(-.39640) 

0.04225 

85.0 

hcoh+h 2 

-1 15. 37823(-. 41465) 

0.03828 

71.1 


a Energy in Eh- The first energy is the CCI energy, while the energy in parenthesis 
includes a multi-reference Davidson correction and is with respect to -115. Eh- 
6 zero-point energy in Eh- 

c relative energy in keal/mol including zero-point energy and a multi-reference 
Davidson’s correction. 
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Table III. Computed saddle point frequencies and r otational constants(cm *)■ 

CH 2 0-H 2 hcoh-h 2 ch 2 .h 2 o ch 2 -h 2 o 


Wi 

3195 

4130 

3860 

3723 

Ul 2 

2295 

3199 

3744 

3319 

0^3 

1740 

2323 

3187 

3222 

0/ 4 

1574 

1564 

3113 

2091 

<^5 

1429 

1498 

1691 

1560 

Ct>6 

1369 

1347 

1515 

1524 

Ct?7 

916 

1291 

802 

1119 

^8 

2877i 

1122 

610 

949 

U>9 

3278 

934 

337 

703 


1273 

622 

334 

438 

u>n 

1211 

513 

161 

387 

t*>12 

1065 

1414t 

26 

1850i 

A 

3.345 

3.042 

4.117 

4.640 

B 

0.944 

0.839 

0.358 

0.496 

C 

0.863 

0.764 

0.349 

0.478 
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Table IV. Computed energies and zero-point corrections' 1 . 


CH 2 + H 2 0 CH 2 - H 2 0 


r co 

CAS Energy 

Cl Energy 6 

zero-point energy 

total c 

3.834 

-115.00924 

-115. 36931 (-.40561) 

0.04366 

(-.36225) 

3.929 

-115.01305 

-115.36889(-.40409) 

0.04449 

(-.35960) 

4.076 

-115.01765 

-115.36968(-.40420) 

0.04550 

(-.35870) 

4.207 

-115.02142 

-115.37055(-.40452) 

0.04546 

(-.35906) 

4.310 

-115.02426 

-1 15. 37130(-. 40480) 

0.04523 

(-.35957) 

4.386 

-115.02635 

- 115. 371 89(-. 40503) 

0.04500 

(-.36003) 

4.437 

-115.02785 

-1 15. 37230(-. 40516) 

0.04482 

(-.36034) 

4.471 

-115.02890 

-115.37254(-.40517) 

0.04470 

(-.36047) 

4.528 

-115.03024 

-115.37237(-. 40450) 

0.04421 

(-.36029) 

5.0 

-115.02970 

-115. 36971 (-.40 136) 



5.5 

-115.02892 

-115. 36727(-. 39843) 



6.0 

-115.02797 

-11 5. 36541 (-.39630) 



6.5 

-115.02710 

-1 15. 36400(-. 39475) 



7.0 

-115.02635 

-1 15. 36294(-. 39358) 



8.0 

-115.02526 

-115. 361 51 (-.39206) 



9.0 

-115.02464 

-115.36076(-.39127) 



10.0 

-115.02433 

-1 15. 36040(-. 39091) 




“ All energies in Eh. r co is in a. u. 
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6 The first energy is the CCI energy, while the energy in parenthesis includes a 
multi-reference Davidson correction and is with respect to -115. Eh- 
c Energy with respect to -115. Eh* The energy includes the zero-point correction 
and a multi-reference Davidson correction. 
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Qxx Qw Q" 


h 2 o 


0.768 -1.987 


1.955 0.032 


CH 2 J A a 0.679 1.187 0.496 -1.683 

- Properties are in a.u. Quadrapole moment is with respect to the center of mass 
The molecule is in the YZ plane with the C 2 axis in the Z direction. 
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Figure Captions. 


Fig. 1. Schematic diagram of the potential energy surface for CH3 + OH. The 
location of the CH 3 + OH asymptote with respect to CH 3 OH is taken from experi- 
ment, while the locations of the CH 2 OH + H aand CH 3 0 + H asymptotes are from 
previous calculations. 

Fig. 2. The potential for 1 CH 2 + H 2 0 from CCI calculations along the CASSCF 
minimum energy path. See the text. 
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